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Abstract 

The peakons discussed here are singular solutions of the dispersionless 
Camassa-Holm (CH) shallow water wave equation in one spatial dimension. 
These are reviewed in the context of asymptotic expansions and Euler- 
Poincare variational principles. The dispersionless CH equation generalizes 
to the EPDiff equation, whose singular solutions are peakon wave fronts 
in higher dimensions. The reduction of these singular solutions of CH and 
EPDiff to canonical Hamiltonian dynamics on lower dimensional sets may 
be understood, by realizing that their solution ansatz is a momentum map, 
and momentum maps are Poisson. 

1 Introduction 

Camassa and Holm [7J discovered the "peakon" solitary traveling wave solution 
for a shallow water wave, 
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whose fluid velocity u is a function of position x on the real line and time t. The 
peakon traveling wave moves at a speed equal to its maximum height, at which 
it has a sharp peak (jump in derivative). Peakons are an emergent phenomenon, 
solving the initial value problem for a partial differential equation derived by an 
asymptotic expansion of Euler's equations using the small parameters of shallow 
water dynamics. Peakons are nonanalytic solitons, which superpose as 



for sets {p} and {q} satisfying canonical Hamiltonian dynamics. Peakons arise 
for shallow water waves in the limit of zero linear dispersion in one dimension. 
Peakons satisfy a partial differential equation (PDE) arising from Hamilton's prin- 
ciple for geodesic motion on the smooth invertible maps (diffeomorphisms) with 
respect to the H 1 Sobolev norm of the fluid velocity. Peakons generalize to higher 
dimensions, as well. We explain how peakons were derived in the context of 
shallow water asymptotics and describe some of their remarkable mathematical 
properties. 

2 Shallow water background for peakons 

Euler's equations for irrotational incompressible ideal fluid motion under gravity 
with a free surface have an asymptotic expansion for shallow water waves that 
contains two small parameters, e and 5 2 , with ordering e > 5 2 . These small 
parameters are e = a /ho (the ratio of wave amplitude to mean depth) and 5 2 = 
{ho/lx) 2 (the squared ratio of mean depth to horizontal length, or wavelength). 
Euler's equations are made non-dimensional by introducing x = l x x' for horizontal 
position, z = h z' for vertical position, t = (l x /c )t' for time, rj — arj' for surface 
elevation and Lp = (gl x a/c )(p' for velocity potential, where c = \/gh is the 
mean wave speed and g is the constant gravity. The quantity a = a' /(hnpcfy is 
the dimensionless Bond number, in which p is the mass density of the fluid and 
a 1 is its surface tension, both of which are taken to be constants. After dropping 
primes, this asymptotic expansion yields the nondimensional Korteweg-de Vries 
(KdV) equation for the horizontal velocity variable u = ip x (x,t) at linear order in 
the small dimensionless ratios e and 5 2 , as the left hand side of 
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Here, partial derivatives are denoted using subscripts, and boundary conditions 
are u = and u x = at spatial infinity on the real line. The famous sech 2 (x — t) 
traveling wave solutions (the solitons) for KdV (|3| arise in a balance between 
its (weakly) nonlinear steepening and its third-order linear dispersion, when the 
quadratic terms in e and 5 2 on its right hand side are neglected. 

In equation ([3]), a normal form transformation due to Kodama [22] has been 
used to remove the other possible quadratic terms of order 0(e 2 ) and 0(S 4 ). The 
remaining quadratic correction terms in the KdV equation ^ may be collected 
at order 0(eS 2 ). These terms may be expressed, after introducing a "momentum 
variable," 

m = u — vb 2 u xx , (4) 
and neglecting terms of cubic order in e and 5 2 , as 

e 5 2 
m t + m x + -(um x + bmu x ) + — (1 - 3a)u xxx = . (5) 
z o 

In the momentum variable m = u — u5 2 u xx , the parameter v is given by [TT] 

19 -30(J -45a 2 

(6) 



60(1 -3a) 



Thus, the effects of 6 2 — dispersion also enter the nonlinear terms. After restor- 
ing dimensions in equation ^ and rescaling velocity u by (6+1), the following 
"6— equation" emerges, 

m t + c m x + um x + b mu x + Tu xxx = , (7) 

where m = u — a 2 u xx is the dimensional momentum variable, and the constants a 2 
and T/co are squares of length scales. When a 2 — > 0, one recovers KdV from the 
b— equation ([7]), up to a rescaling of velocity. Any value of the parameter 6^—1 
may be achieved in equation (J7|) by an appropriate Kodama transformation [TT] . 

As we have emphasized, the values of the coefficients in the asymptotic analysis 
of shallow water waves at quadratic order in their two small parameters only hold, 
modulo the Kodama normal-form transformations. Hence, these transformations 
may be used to advance the analysis and thereby gain insight, by optimizing the 
choices of these coefficients. The freedom introduced by the Kodama transfor- 
mations among asymptotically equivalent equations at quadratic order in e and 
S 2 also helps to answer the perennial question, "Why are integrable equations so 
ubiquitous when one uses asymptotics in modeling?" 
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Integrable cases of the 6— equation ([7]). The cases 6 = 2 and 6 = 3 are special 
values, for which the 6— equation becomes a completely integrable Hamiltonian 
system. For 6 = 2, equation §f§ specializes to the integrable CH equation of 
Camassa and Holm [TJ. The case 6 = 3 in (|7|) recovers the integrable DP equation 
of Degasperis and Procesi [UJ . These two cases exhaust the integrable candidates 
for ([7]), as was shown using Painleve analysis. The 6— family of equations ([7]) was 
also shown in [26] to admit the symmetry conditions necessary for integrability, 
only in the cases 6 = 2 for CH and 6 = 3 for DP. 

The 6— equation ^ with 6 = 2 was first derived in Camassa and Holm [7J 
by using asymptotic expansions directly in the Hamiltonian for Euler's equations 
governing inviscid incompressible flow in the shallow water regime. In this anal- 
ysis, the CH equation was shown to be bi-Hamiltonian and thereby was found 
to be completely integrable by the inverse scattering transform (1ST) on the real 
line. Reviews of 1ST may be found, for example, in Ablowitz et al. [JJ, Dubrovin 
et al. [10], Novikov et al. [27]. For discussions of other related bi-Hamiltonian 
equations, see [9J. 

Camassa and Holm [7J also discovered the remarkable peaked soliton (peakon) 
solutions of ( 

6 = 2. The peakons arise as solutions of ([7]), when Co = and T = in the absence 
of linear dispersion. Peakons move at a speed equal to their maximum height, at 
which they have a sharp peak (jump in derivative). Unlike the KdV soliton, the 
peakon speed is independent of its width (a). Periodic peakon solutions of CH 
were treated in Alber et al. [2]. There, the sharp peaks of periodic peakons were 
associated with billiards reflecting at the boundary of an elliptical domain. These 
billiard solutions for the periodic peakons arise from geodesic motion on a tri-axial 
ellipsoid, in the limit that one of its axes shrinks to zero length. 

Before Camassa and Holm derived their shallow water equation in [7J , a class 
of integrable equations existed, which was later found to contain equation ([7]) with 
6 = 2. This class of integrable equations was derived using hereditary symmetries 
in Fokas and Fuchssteiner [13]. However, equation (|7| was not written explicitly, 
nor was it derived physically as a shallow water equation and its solution properties 
for 6 = 2 were not studied before Camassa and Holm [7J. See Fuchssteiner [T5] 
for an insightful history of how the shallow water equation ^ in the integrable 
case with 6 = 2 relates to the mathematical theory of hereditary symmetries. 

Equation ([7]) with 6 = 2 was recently re-derived as a shallow water equation 
by using asymptotic methods in three different approaches in Dullin et al. [TT] . 
in Fokas and Liu [H] and also in Johnson [2JJ . These three derivations all used 
different variants of the method of asymptotic expansions for shallow water waves 



lp) for the CH equation on the real line, given by (TH) in the case 
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in the absence of surface tension. Only the derivation in Dullin et al. [TT] used the 
Kodama normal-form transformations to take advantage of the non-uniqueness of 
the asymptotic expansion results at quadratic order. 

The effects of the parameter b on the solutions of equation (j7| were investigated 
in Holm and Staley [19], where b was treated as a bifurcation parameter, in the 
limiting case when the linear dispersion coefficients are set to Cq = and T = 0. 
This limiting case allows several special solutions, including the peakons, in which 
the two nonlinear terms in equation ([7]) balance each other in the absence of linear 
dispersion. 

3 Peakons: Singular solutions without linear 
dispersion in one spatial dimension 

Peakons were first found as singular soliton solutions of the completely integrable 
CH equation. This is equation ^ with 6 = 2, now rewritten in terms of the 
velocity, as 

v* + CqU x + 3uu x + Tu xxx 

= a 2 [u xxt + 2u x u xx + uu xxx ) . (8) 

Peakons were found in [7j to arise in the absence of linear dispersion. That is, they 
arise when Co = and T = in CH ([8]). Specifically, peakons are the individual 
terms in the peaked iV— soliton solution of CH (|8j) for its velocity, 

N 

u (x, t) = J2 Pb{t)e~ lx - qb{t)l/a , (9) 

6=1 

in the absence of linear dispersion. Each term in the sum is a solition with a 
sharp peak at its maximum. Hence, the name "peakon." Expressed using its 
momentum, m = (1 — a 2 d x )u, the peakon velocity solution ^ of dispersionless 
CH becomes a sum over a delta functions, supported on a set of points moving 
on the real line. Namely, the peakon velocity solution ^ implies 

N 

m(x,t) = 2a Pb(t)S(x - q b (t)) , (10) 

6=1 

because of the relation (1 — a 2 dl)e~^ x ^ a = 2a5(x). These solutions satisfy the 
b— equation Q for any value of 6, provided Co = and T = 0. 
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Thus, peakons are singular momentum solutions of the dispersionless 
b— equation, although they are not stable for every value of b. From numerical 
simulations [19], peakons are conjecture to be stable for b > 1. In the integrable 
cases b = 2 for CH and b = 3 for DP, peakons are stable singular soliton solu- 
tions. The spatial velocity profile e~^ x ^ a /(2a) of each separate peakon in (JoJ) is 
the Green's function for the Helmholtz operator on the real line, with vanishing 
boundary conditions at spatial infinity. Unlike the KdV soliton, whose speed and 
width are related, the width of the peakon profile is set by its Green's function, 
independently of its speed. 



Integrable peakon dynamics of CH. Substituting the peakon solution ansatz 



9|) and (10) into the dispersionless CH equation, 
m t + um x + 2mu x = , with 



m 



u 



a 2 u,, 



(11) 



yields Hamilton's canonical equations for the dynamics of the discrete set of peakon 
parameters p a (t) and q a (t), 



. /,x dh N . dh N 

q a {t) = — and p a (t) = - -r— 
dp a dq a 

for a = 1, 2, . . . , N, with Hamiltonian given by [7j, 

N 



(12) 



>N 



\ paPb ' 



-\<la-qb\/ai 



(13) 



a,b=l 



Thus, one finds that the points x = q a (t) in the peakon solution ^ move with 
the flow of the fluid velocity u at those points, since u(q a (t),t) = q a (t). This 
means the q a (t) are Lagrangian coordinates. Moreover, the singular momentum 



solution ansatz (10) is the Lagrange-to-Euler map for an invariant manifold of the 



dispersionless CH equation ( 11 ). On this finite-dimensional invariant manifold for 



the partial differential equation (11), the dynamics is canonically Hamiltonian. 



With Hamiltonian (13), the canonical equations (12) for the 2iV canonically 



conjugate peakon parameters p a (t) and q a (f) were interpreted in [7j as describing 
geodesic motion on the N— dimensional Riemannian manifold whose co-metric is 



9 tJ ({q}) = e 



-\<n-ij\/ a 



Moreover, the canonical geodesic equations arising from 



Hamiltonian (13) comprise an integrable system for any number of peakons N. 



This integrable system was studied in [7] for solutions on the real line, and in 
[21 125] and references therein, for spatially periodic solutions. 
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Being a completely integrable Hamiltonian soliton equation, the continuum 
CH equation (l8j) has an associated isospectral eigenvalue problem, discovered in 
[7] for any values of its dispersion parameters Co and T. Remarkably, when cq = 
and T = 0, this isospectral eigenvalue problem has a purely discrete spectrum. 
Moreover, in this case, each discrete eigenvalue corresponds precisely to the time- 
asymptotic velocity of a peakon. This discreteness of the CH isospectrum in the 
absence of linear dispersion implies that only the singular peakon solutions (10) 
emerge asymptotically in time, in the solution of the initial value problem for the 
dispersionless CH equation (11). This is borne out in numerical simulations of 
the dispersionless CH equation (11), starting from a smooth initial distribution of 
velocity [III US- 
Figure [T] shows the emergence of peakons from an initially Gaussian velocity 
distribution and their subsequent elastic collisions in a periodic one-dimensional 
domain]]] This figure demonstrates that singular solutions dominate the initial 
value problem and, thus, that it is imperative to go beyond smooth solutions for 
the CH equation; the situation is similar for the EPDiff equation. 



Peakons as mechanical systems. Being governed by canonical Hamiltonian 
equations, each N— peakon solution can be associated with a mechanical system 
of moving particles. Calogero et al. [6j further extended the class of mechan- 
ical systems of this type. The r-matrix approach was applied to the Lax pair 
formulation of the N— peakon system for CH by Ragnisco and Bruschi [28], who 
also pointed out the connection of this system with the classical Toda lattice. 
A discrete version of the Adler-Kostant-Symes factorization method was used by 
Suris [29] to study a discretization of the peakon lattice, realized as a discrete 
integrable system on a certain Poisson submanifold of gl(N) equipped with an 
r-matrix Poisson bracket. Beals et al. |5J used the Stieltjes theorem on continued 
fractions and the classical moment problem for studying multi-peakon solutions 
of the CH equation. Generalized peakon systems are described for any simple Lie 
algebra by Alber et al. [2]. 



Pulsons: Generalizing the peakon solutions of the dispersionless 
b— equation for other Green's functions. The Hamiltonian in equation 
(13) depends on the Green's function for the relation between velocity u and mo- 
mentum m. However, the singular momentum solution ansatz (10) is independent 
of this Green's function. Thus, as discovered in Fringer and Holm |16j . 

lr rhe figures in this article were kindly supplied by Martin Staley 
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Figure 1: A smooth localized (Gaussian) initial condition for the CH equation 
breaks up into an ordered train of peakons as time evolves (the time direction 
being vertical). The peakon train eventually wraps around the periodic domain, 
thereby allowing the leading peakons to overtake the slower emergent peakons 
from behind in collisions that cause phase shifts, as discussed in [TJ. 



The singular momentum solution ansatz (10) for the dispersionless equation, 

m t + um x + 2mu x = , with u = g * m , (14) 

provides an invariant manifold on which canonical Hamiltonian dynamics occurs, 
for any choice of the Green's function g relating velocity u and momentum m by 
the convolution u = g * m. 

The fluid velocity solutions corresponding to the singular momentum ansatz 
(10) for equation (14) are the pulsons. Pulsons are given by the sum over N 
velocity profiles determined by the Green's function g, as 



N 



u(x,t) = ^2 Pa(t)g(x,q a (t)) 



(15) 



a=l 



Again for (14), the singular momentum ansatz (10) results in a finite-dimensional 
invariant manifold of solutions, whose dynamics is canonically Hamiltonian. The 
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Hamiltonian for the canonical dynamics of the 27V" parameters p a (t) and q a (t) in 



the "pulson" solutions ( 15 ) of equation ( 14 ) is 



1 N 



(16) 



a,b=l 



Again for the pulsons, the canonical equations for the invariant manifold of sin- 
gular momentum solutions provide a phase-space description of geodesic motion, 
this time with respect to the co- metric given by the Green's function g. Mathe- 
matical analysis and numerical results for the dynamics of these pulson solutions 
are given in [16J. These results describe how the collisions of pulsons (15) depend 
upon their shape. 



Compactons in the I /a 2 — > limit of CH. As mentioned earlier, in the limit 
that a 2 — > 0, the CH equation ^ becomes the KdV equation. In the opposite 
limit that I /a 2 — > CH becomes the Hunter-Zheng equation [20] 

(ut + uu x ) xx = ^{ul) x (Hunter-Zheng) 

This equation has "compacton" solutions, whose collision dynamics was studied 
numerically and put into the present context in [16]. The corresponding Green's 
function satisfies —d 2 g(x) = 25(x), so it has the triangular shape, g(x) — 1 — \x\ 
for | a; | < 1, and vanishes otherwise, for \x\ > 1. That is, the Green's function in 
this case has compact support; hence, the name "compactons" for these pulson 
solutions, which as a limit of the integrable CH equations are true solitons, solvable 
by 1ST. 



Pulson solutions of the dispersionless b— equation. Holm and Staley [19] 
give the pulson solutions of the traveling wave problem and their elastic collision 
properties for the dispersionless b— equation, 

m t + um x + b mu x = , with u = g * m , (17) 

with any (symmetric) Green's function g and for any value of the parameter b. 
Numerically, pulsons and peakons are both found to be stable for b > 1, [19]. The 
reduction to noncanonical Hamiltonian dynamics for the invariant manifold of 



singular momentum solutions ( 10 ) of the other integrable case 6 = 3 with peakon 
Green's function g(x,y) = e~' a; ^' a j s found in [9]. 
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4 Euler-Poincare theory in more dimensions 



Generalizing the peakon solutions of the CH equation to higher di- 
mensions. In [19J, weakly nonlinear analysis and the assumption of columnar 
motion in the variational principle for Euler's equations are found to produce the 



two-dimensional generalization of the dispersionless CH equation (11). This gen- 



eralization is the Euler-Poincare (EP) equation [18] for the Lagrangian consisting 
of the kinetic energy, 



|u| 2 + a 2 (divu) dxdy . 



in which the fluid velocity u is a two-dimensional vector. Evolution generated 
by kinetic energy in Hamilton's principle results in geodesic motion, with respect 
to the velocity norm ||u||, which is provided by the kinetic energy Lagrangian. 
For ideal incompressible fluids governed by Euler's equations, the importance of 
geodesic flow was recognized by Arnold [3] for the L 2 norm of the fluid velocity. 
The EP equation generated by any choice of kinetic energy norm without imposing 
incompressibility is called "EPDiff," for "Euler-Poincare equation for geodesic 
motion on the diffeomorphisms." EPDiff is given by 



+ u • V j m + Vu T • m + m(div u) 







(19) 



with momentum density m = 5£/5u, where £ = |||u|| 2 is given by the kinetic 
energy, which defines a norm in the fluid velocity ||u||, yet to be determined. By 
design, this equation has no contribution from either potential energy, or pressure. 
It conserves the velocity norm ||u|| given by the kinetic energy. Its evolution 
describes geodesic motion on the diffeomorphisms with respect to this norm 



An alternative way of writing the EPDiff equation (19) in either two, or three 
dimensions is, 



0_ 
Of 



m — u x curl m + V(u • m) + m(div u) = . 



(20) 



This form of EPDiff involves all three differential operators, curl, gradient and 
divergence. For the kinetic energy Lagrangian £ given in (18), which is a norm 
for irrotational flow (with curlu = 0), we have the EPDiff equation (19) with 
momentum m = 5£/5u = u — « 2 V(divu). 



EPDiff (19) may also be written intrinsically as 
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where ad* is the L 2 dual of the ad-operation (commutator) for vector fields. See jH 
EI] for additional discussions of the beautiful geometry underlying this equation. 



Reduction to the dispersionless CH equation in ID. In one dimension, 



the EPDiff equation (19 21) with Lagrangian i given in (18) simplifies to the 



dispersionless CH equation (11). The dispersionless limit of the CH equation 



appears, because we have ignored potential energy and pressure. 



Strengthening the kinetic energy norm to allow for circulation. The 



kinetic energy Lagrangian (18) is a norm for irrotational flow, with curlu 



0. However, inclusion of rotational flow requires the kinetic energy norm to be 
strengthened to the H\ norm of the velocity, defined as 



|u| 2 + a 2 (divu) 2 + a 2 (curl u)' 

1„ „, 



dxdy 



lur + a Vu 



dxdy 



u 



(22) 



Here we assume boundary conditions that give no contributions upon integrating 
by parts. The corresponding EPDiff equation is (19) with m = 5£/5u = u — 



a Au . This expression involves inversion of the familiar Helmholtz operator in the 
(nonlocal) relation between fluid velocity and momentum density. The H\ norm 
|| u ||^i for the kinetic energy (22) also arises in three dimensions for turbulence 
modeling based on Lagrangian averaging and using Taylor's hypothesis that the 
turbulent fluctuations are "frozen" into the Lagrangian mean flow [12]. 



Generalizing the CH peakon solutions to n dimensions. Building on the 
peakon solutions ^ for the CH equation and the pulsons (JTsJ) for its generalization 
to other traveling- wave shapes in jTB], Holm and Staley [TJ5] introduced the fol- 
lowing measure- valued singular momentum solution ansatz for the n— dimensional 



solutions of the EPDiff equation (19): 



N 

m(x, t) = 22 pa ( s ' *) 5 ( x ~ Qa ( s > *) ) ds - 

a=l ^ 



(23) 



These singular momentum solutions, called "diffeons," are vector density functions 
supported in M n on a set of iV surfaces (or curves) of codimension (n — k) for 
sei 4 with k < n. They may, for example, be supported on sets of points (vector 
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peakons, k = 0), one- dimensional filaments (strings, k = 1), or two-dimensional 
surfaces (sheets, = 2) in three dimensions. 

Figure [2] shows the results for the EPDiff equation when a straight peakon 
segment of finite length is created initially moving rightward (East). Because of 
propagation along the segment in adjusting to the condition of zero speed at its 
ends and finite speed in its interior, the initially straight segment expands outward 
as it propagates and curves into a peakon "bubble." 




Figure 2: A peakon segment of finite length is initially moving rightward (East). 
Because its speed vanishes at its ends and it has fully two-dimensional spatial 
dependence, it expands into a peakon "bubble" as it propagates. (The colors 
indicate speed: red is highest, yellow is less, blue low, grey zero.) 

Figure [3] shows an initially straight segment whose velocity distribution is 
exponential in the transverse direction, but is wider than a for the peakon solution. 
This initial velocity distribution evolves under EPDiff to separate into a train of 
curved peakon "bubbles," each of width a. This example illustrates the emergent 
property of the peakon solutions in two dimensions. This phenomenon is observed 
in nature, for example, as trains of internal wave fronts in the south China Sea 
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Figure 3: An initially straight segment of velocity distribution whose exponential 
profile is wider than the width a for the peakon solution will break up into a 
train of curved peakon "bubbles," each of width a. This example illustrates the 
emergent property of the peakon solutions in two dimensions. 



Substitution of the singular momentum solution ansatz (23) into the EPDiff 



equation (19) implies the following integro-partial-differential equations (IPDEs) 
for the evolution of the parameters {P} and {Q}, 



0_ 
Of 

d_ 
Of 



N 

Q"M) = E / P\s',t)G(Q a (s,t)-Q\s',t))ds'. 

6=1 

N 

P a (s,t) = -J2 / (P a (s,t).P b (s',t)) 

6=1 J 

G(Q a (s,t) -Q b (s',t))ds' . 



(24) 



dQ a (s,t) 



Importantly for the interpretation of these solutions, the coordinates s G M. h 
turn out to be Lagrangian coordinates. The velocity field corresponding to the 
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Figure 4: A single collision is shown involving reconnection as the faster peakon 
segment initially moving Southeast along the diagonal expands, curves and 
obliquely overtakes the slower peakon segment initially moving rightward (East). 
This reconnection illustrates one of the collision rules for the strongly two- 
dimensional EPDiff flow. 



momentum solution ansatz (23) is given by 



u(x, t) = G * m 



N 

£ / P b (s',t)G(x-Q h (s',t))ds' 

6=1 J 



(25) 



for u G W 1 . When evaluated along the curve x = Q a (s,t), this velocity satisfies, 



N 

u(Q a (s,t),t) = / P b (s',t)G(Q a (s,t)-Q b (s',t))ds' 

6=1 •' 



dQ a (s,t) 
dt 



(26) 
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Consequently, the lower-dimensional support sets defined on x = Q a (s,t) and 

i fc move with the fluid velocity. This means 



parameterized by coordinates s G 
the sGl' are Lagrangian coordinates. Moreover, equations (24) for the evolution 
of these support sets are canonical Hamiltonian equations, 



dt 4 1 ' } 5P- 



0_ 
Of 



P a (s,t) 



The corresponding Hamiltonian function H 



N 



6H 



N 



5Q a 



x R n ) N -> R is, 



(27) 



hn = Iff it, (^ a (s,t)-p h (s',t)) 



a , b=l 



G(Q a {s,t),Qj{s',t))dsds' . 



(2f 



This is the Hamiltonian for geodesic motion on the cotangent bundle of a set of 
curves Q a (s, t) with respect to the metric given by G. This dynamics was investi- 
gated numerically in [19] to which we refer for more details of the solution prop- 
erties. One important result found numerically in [TH] is that only codimension- 
one singular momentum solutions appear to be stable under the evolution of the 
EPDiff equation. Thus, we have 

Stability for codimension-one: the singular momentum solutions of EPDiff are 
stable, as points on the line (peakons), as curves in the plane (filaments, or wave 
fronts), or as surfaces in space (sheets). 

Proving this stability result analytically remains an outstanding problem. The 
stability of peakons on the real line is proven in [Hj. 



Reconnections in oblique overtaking collisions of peakon wave fronts. 

Figures [4] and [5] show results of oblique wave front collisions producing reconnec- 
tions for the EPDiff equation in two dimensions. Figure [4] shows a single oblique 
overtaking collision, as a faster expanding peakon wave front overtakes a slower 
one and reconnects with it at the collision point. Figure [5] shows a series of recon- 
nections involving the oblique overtaking collisions of two trains of curved peakon 
filaments, or wave fronts. 



The peakon reduction is a momentum map. As shown in [T7], the singular 
solution ansatz (23 ) is a momentum map from the cotangent bundle of the smooth 
embeddings of lower dimensional sets W G M. n , to the dual of the Lie algebra of 
vector fields defined on these sets. (Momentum maps for Hamiltonian dynamics 
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Figure 5: A series of multiple collisions is shown involving reconnections as the 
faster wider peakon segment initially moving Northeast along the diagonal ex- 
pands, breaks up into a wave train of peakons, each of which propagates, curves 
and obliquely overtakes the slower wide peakon segment initially moving right- 
ward (East), which is also breaking up into a train of wavefronts. In this series 
of oblique collision, the now-curved peakon filaments exchange momentum and 
reconnect several times. 



are reviewed in [24], for example.) This geometric feature underlies the remark- 
able reduction properties of the EPDiff equation, and it also explains why the 
reduced equations must be Hamiltonian on the invariant manifolds of the singular 
solutions; namely, because momentum maps are Poisson maps. This geometric 



feature also underlies the singular momentum solution (23) and its associated ve- 



locity ( 25 ) which generalize the peakon solutions, both to higher dimensions and 



to arbitrary kinetic energy metrics. The result that the singular solution ansatz 



(23) is a momentum map helps to organize the theory, to explain previous results 



and to suggest new avenues of exploration. 
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